Unexpected spatial distribution of bubble rearrangements in coarsening foams 
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Foams are ideal model systems to study stress-driven dynamics, as stress-imbalances within the 
system are continuously generated by the coarsening process, which unlike thermal fluctuations, can 
be conveniently quantified by optical means. However, the high turbidity of foams generally hinders 
the detailed study of the temporal and spatial distribution of rearrangement events, such that definite 
assessments regarding their contribution to the overall dynamics could not be made so far. In this 
paper, we use novel light scattering techniques to measure the frequency and position of events within 
a large sample volume. As recently reported (A. S. Gittings and D. J. Durian, Phys. Rev. E, 2008, 
78, 066313), we find that the foam dynamics is determined by two distinct processes: intermittent 
bubble rearrangements of finite duration and a spatially homogeneous quasicontinuous process. Our 
experiments show that the convolution of these two processes determines the age-dependence of the 
mean dynamics, such that relations between intermittent rearrangements and coarsening process 
can not be established by considering means. By contrast the use of the recently introduced photon 
correlation imaging technique (A. Duri, D. A. Sessoms, V. Trappe, and L. Cipelletti, Phys. Rev. 
Lett., 2009, 102, 085702) enables us to assess that the event frequency is directly determined by 
the strain-rate imposed by the coarsening process. Surprisingly, we also find that, although the 
distribution of successive events in time is consistent with a random process, the spatial distribution 
of successive events is not random: rearrangements are more likely to occur within a recently 
rearranged zone. This implies that a topological rearrangement is likely to lead to an unstable 
configuration, such that a small amount of coarsening-induced strain is sufficient to trigger another 
event . 

PACS numbers: 47.60. Dx 47.55.D- 47.20.Ky 



INTRODUCTION 



Foams consist of gas bubbles that are densely packed 
in a continuous liquid phase. Because of the high packing 
fraction, the bubbles are deformed, such that the system 
is constantly under internal stress. Though thermal mo- 
tion is insufficient to lead to any significant restructuring 
of the random bubble configuration, foams reconfigure in 
time because of coarsening: larger bubbles growing at the 
expense of the smaller ones due to the difference in their 
Laplace pressures. This process leads to imbalances of 
the internal stresses, which are released by intermittent, 
local rearrangements of bubbles [HQ. Thus, the coars- 
ening process provides the system with an internal me- 
chanical means of thermalization that eventually allows 
for a complete reconfiguration of the foam. Such stress- 
driven relaxation is the hallmark of a variety of other soft 
glassy systems [U , such as colloidal gels Q , concentrated 
microgel systems and surfactant pha7es@; however, 
in these systems, the mechanism leading to imbalanccd 
stresses is not well understood. By contrast, in foams 
the coarsening process is a quantifiable measure of the 
temporal evolution of stress-imbalances, and thus foams 
are ideal benchmark systems to study stress driven re- 
laxation processes. 
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Their dynamics has been previously investigated using 
diffusing- wave spectroscopy (DWS) @, 043], an exten- 
sion of dynamic light scattering to the multiple scattering 
limit 13, [l4|. However, traditional DWS relies on exten- 
sive time and space averaging, thereby precluding direct 
characterization of the intermittent dynamics. Indeed, 
recent speckle visibility experiments revealed that the dy- 
namics of foams exhibits features that were not consid- 
ered in previous modelling of the time-averaged correla- 
tion functions (l5|. In particular, it was shown that the 
dynamic light scattering signal was not only determined 
by intermittent rearrangements of the local bubble con- 
figuration, but that another process of undefined origin 
was contributing as well. Such mix of dynamic processes, 
however, significantly complicates the interpretation of 
time and space averaged dynamic light scattering signals 
that have been used so far to characterize foam dynam- 
ics. In particular, it obscures the evaluation of the re- 
lation between the intermittent restructuring events and 
the coarsening process, which is essential to understand 
the interplay between stress accumulation, local yielding 
and structural relaxation. 

In this article, we extend the investigations of foam dy- 
namics to spatially resolved dynamics, applying the re- 
cently introduced photon correlation imaging scheme [l6| 
to the backscattering plane of a coarsening foam. To cor- 
relate our findings to previous work, we simultaneously 
perform time resolved experiments [T3, EH , in the classi- 
cal transmission far field geometry of DWS. Our results 
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reveal that the event frequency is directly correlated to 
the strain imposed by the coarsening process, which is not 
evident when considering the time and spatially averaged 
signals within the frame of previous models, where the 
existence of an additional dynamic process was not con- 
sidered. Moreover, we find that though the occurrence of 
events is random in time it is not in space. Indeed, our 
results suggest that successive events preferentially occur 
at the same location. 



II. EXPERIMENTAL 

Our foam is a commercial Gillette shaving foam that 
consists of polydispersc bubbles tightly packed in an 
aqueous solution of stearic acid and triethanolamine. 
As demonstrated in a number of investigations 0, 0- 
03 El EE El-il, this type of foam exhibits dynami- 
cal and rheological properties that are remarkably insen- 
sitive to variations in the detailed composition, which 
inevitably will vary as a function of the market evolu- 
tion. This insensitivity denotes Gillette shaving foams as 
an ideal benchmark system to study different features of 
foam dynamics, yielding results that can be compared to 
the large number of prior investigations. 

To characterize the coarsening process of our system, 
we perform both optical microscopy and diffuse light 
transmittance experiments fl3j . Immediately after pro- 
duction, the foam is injected into a rectangular cell of 
thickness L = 2mm or 5mm and kept at a stable ambi- 
ent temperature of T = 21°C. The age of the sample t w 
is defined as the time elapsed since sample production. 
We record the evolution of the bubble size by imaging 
the foam surface, where we determine the mean bub- 
ble size (at m } based on the number of bubbles detected 
within the area defined by our field of view. As previ- 
ously shown, (a tul ) is directly proportional to the trans- 
port mean free path I* that we measure in diffuse trans- 
mittance experiments 

HEH. To perform these experi- 
ments we illuminate our sample with a collimated laser 
beam with a wavelength of 532 nm in vacuo and a di- 
ameter of 9 mm. A schematic layout of our experimental 
set-up is shown in Fig. Q] We detect the mean inten- 
sity of the transmitted light as a function of t w using a 
photomultiplier and determine the age dependent trans- 
port mean free path of our sample using the transmitted 
intensity of a reference of known I* to calibrate our exper- 
iment. These experiments are performed simultaneously 
to the diffusing wave spectroscopy experiments outlined 
below, which enables us to determine slight variations 
in the coarsening process in different runs of the experi- 
ment. By mapping the evolution of /* to the evolution of 
(at w ), we find that for our foam I*/ (2 • a tw ) = 4.6 ± 0.3, 
which is somewhat larger than the typical values found 
in previous experiments probing the optical properties of 
Gillette foams 0, Ell ■ In agreement with previous results 
we find that the bubble size evolution is well described 
by (at m ) 2 — (ao) 2 oc /J 1 2 — /q 2 cx t w , with (ao) and Iq 
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FIG. 1: Schematic of light scattering set-up used to record si- 
multaneously the speckle pattern formed in the far-field trans- 
mission geometry and in the backscattering plane of the sam- 
ple using CCD-cameras. To assess the evolution of the bubble 
size at each run, the transmitted light intensity is measured 
continuously using a photo-multiplier (PM) detection scheme. 



respectively the bubble radius and the transport mean 
free path at t w = 0, ■ Our main experiments are 
performed at two different ages, min (age 1) 

and t w = 663 min (age 2), where t w denotes the age of 
the sample at the middle of the experiment. The bubble 
radii are respectively ~ 50/im and ~ 76/im at age 1 and 
age 2. The duration of our experiments is 800s at age 1 
and 1200s at age 2, which ensures that the bubble size 
does not change by more than 2% in the course of our 
experiments. 

To characterize the dynamics of our foam at age 1 and 
age 2, we take advantage of the concepts of diffusing wave 
spectroscopy. Of particular relevance is here the model 
that has been developed to describe a gradual reconfig- 
uration of a foam by intermittent rearrangements of the 
local bubble configuration 0, [l!| . In this model, the re- 
arrangements are presumed to occur randomly in space 
and time, to be instantaneous and to affect a well defined 
localized region of radius £. Assuming that all scattered 
light passing a rearranged region is completely dephased, 
while photon paths not passing through a rearranged re- 
gion remain unchanged, the field autocorrelation gi (r) 
can be regarded as a direct measure of the fraction of 
paths not rearranged during the lag time r; for transmis- 
sion experiments gi (r) is then well approximated by 

gl (T)^cxp{-R-lne-(L/n 2 -r} (1) 

with R the rate of intermittent rearrangements per unit 
volume. The characteristic decay time of the field corre- 
lation function (i? • §tt£ 3 ) thus denotes the mean time 
between events at a given location. To gain a physical 
intuition about the impact of a single event on gi (r) for 
a given transmission experiment, we can rewrite equ. [1] 
with respect to the actual volume element observed, the 
scattering volume V sc . 



H (t) = exp <^ -r 



v 3a 



Vsc 



(L/rf 



(2) 



3 



Here, Ty sc is the rate of events occurring within the scat- 
tering volume, such that Ty 3r ■ r corresponds to the num- 
ber of events occurring within V sc during r. The fraction 
of the paths randomized per event is (^ir£ 3 /V sc ) -(L/l*) 2 ; 
it depends on the fraction of the scattering volume rear- 
ranged per event, §tt£ 3 , anci , due to the effective diffusion 
of the light through the sample, on (L/l*) 2 . 

Our experiments are based on the classical schemes of 
dynamic light scattering and on the use of charge coupled 
device (CCD) cameras to capture the intensity fluctua- 
tions of a large number of independent speckles simulta- 
neously. We record the evolution of the speckle pattern 
in time using two different schemes. 

Placing a CCD-detector far from the light-exit-planc 
of our sample, as shown in Fig. [TJ we record the speckle 
pattern in the far-field transmission geometry and pro- 
cess the images using the time resolved correlation (TRC) 
scheme [13, Ell- Here, we calculate a time resolved degree 
of correlation between two speckle patterns recorded at 
time t w and t w + r, according to 



c/ (t w ,r) 



1 



(I p {t w )I p (t w +t)) 



/3 \(Ip(t w )) p (Ip(t w +r)) p 



(3) 



where I p is the intensity recorded at the p-th. pixel of 
the CCD array and (. . .) denotes an average over all 
pixels of the speckle image; (3 is a prefactor that de- 
pends on the speckle-to-pixel size ratio, chosen so that 
c i (twi T ) — 1 f° r t — >• 0, where 7TT denotes a time av- 
erage over the duration of the experiments. Note that 
the time averaged signal corresponds to the usual inten- 
sity correlation function g 2 (r) — 1 = ci(t w ,r), which 
relates to the field correlation function via the Siegert 
relation, g 2 (r) — 1 = gi (r) 2 . The TRC-method enables 
us to access the temporal fluctuations in the dynamics of 
our samples with a spatial average taken over the entire 
scattering volume. The magnitude of the fluctuations 
due to the occurrence of local rearrangement events will 
here depend on the effectiveness of a single event to con- 
tribute to a decrease in cj (t w , r) , which we expect to 
depend on the fraction of the paths dephased per event; 
based on the model discussed above this corresponds to 

^e/Vsc (L/n 2 . 

The second method used is the recently introduced 
photon correlation imaging technique [16j, which is 
closely related to other speckle imaging and near field 
scattering techniques |23l - t2a |. It enables us to resolve 
dynamics with both temporal and spatial resolution. In- 
stead of recording the speckle pattern in the far field, 
where each speckle contains phase contributions of all 
scatterers within the sample, we image the scattered in- 
tensity pattern at the entrance plane of the sample onto 
the camera, as shown in Fig. [TJ and calculate the instan- 
taneous degree of correlation with spatial resolution 
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FIG. 2: Representative dynamic activity maps depicting re- 
spectively a moment at which an event occurs [(a) age 1 and 
(c) age 2] and a moment during the periods between events 
[(b) age 1 and (d) age 2] . The scale bar corresponds to 1 mm. 
Note that the spatial extent of the rearrangement typically 
scales with the bubble size. 



where (. . .) ROI denotes an average over pixels belong- 
ing to a region of interest (ROI) centred around the po- 
sition s. The spatially resolved dynamics of the foam 
at a moment defined by t w and t w + t can then be 
presented graphically by constructing dynamical activ- 
ity maps, where each pixel of the map is colour coded to 
represent the local degree of correlation cf' (t w ,r) char- 
acterizing the fluctuations of diffuse light emanating from 
a specific area of the sample. Examples of such dynami- 
cal activity maps are shown in Fig. [3J The occurrence of 
events is clearly depicted as a bright patch correspond- 
ing to a low Cj (t w , t) in the dynamical activity map, as 
shown for age 1 and age 2 in Fig. [5fa) and (c), respec- 
tively. As our samples are turbid, the appearance of a 
spatially resolved event will be limited to some critical 
depth within the sample, since the backscattered light is 
effectively diffusing on its way back to the entrance plane 
of the cell. To address this issue, we conduct combined 
backscattering and transmission experiments where we 
image respectively the scattered intensity pattern at the 
entrance and exit plane of the cell onto two identical cam- 
eras, aligning the optics so that the observation windows 
are the same in both geometries. We presume that spa- 
tial resolution of an event in transmission can only be 
obtained when the event occurs near the exit plane of the 
cell. Thus, if for a given cell thickness, a spatially distinct 
event in transmission can also be resolved in backscatter- 
ing, we assess that the critical depth for the visibility of 
spatially resolved events is beyond the thickness L. Con- 
ducting different experiments varying L, we assess that 
the critical depth corresponds to L/l* ~ 7.5, which is 
broadly consistent with estimations probing the resolu- 
tion depth in speckle imaging techniques [27]. 

Clearly, in both our experiments, the transmission far 
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field and the backscattcring photon correlation imaging 
experiment, the visibility of events depends on the thick- 
ness of the cell. To obtain the maximal visibility of the 
events, while still maintaining L/l* large enough to be 
within the limit of the multiple scattering approxima- 
tions [lij , we perform our experiments at age 1 in a cell 
with L = 2 mm and at age 2 in a cell with L = 5 mm. 
With I* = 0.44 mm for age 1 and I* = 0.73 mm for age 2 
this corresponds to L/l* = 4.5 (age 1) and L/l* = 6.8 
(age 2). With these L/l* values, we ensure that all 
events occurring throughout the depth of the sample are 
detected in our backscattcring photon correlation imag- 
ing experiments and that the effect of a single event on 
ci {t w , r) measured in our far field transmission TRC- 
experiments is maintained within the same range at both 
ages; here we recall that the effect of a single event on 
Ci (t w ,r) scales with (t; 3 /V sc ) (L/l*) 2 , where V sc can be 
modelled by a truncated cone [13[ and £ scales with the 
bubble size. To be able to directly compare the tempo- 
ral evolution of our signals, we perform both experiments 
simultaneously; in addition, the diffuse transmittancc is 
measured continuously from t w — s to t w = 40000 s, 
which enables us to assess the coarsening process for each 
run. Our CCD detectors comprise 633 x 483 pixels and 
our linear speckle size is 3 pixels in the transmission and 1 
pixel in the backscattering experiments. Before each ex- 
periment we determine the exact magnification M of the 
image taken in the backscattering experiment by imaging 
a finely marked grid that is placed at the entrance plane 
of our cell. Our imaging lens is placed so that M = 1, 
such that the observation area of our experiment is ap- 
proximately 6.3 x 4.8mm 2 . Linear polarizers are placed 
in front of both cameras so that only light with a po- 
larization perpendicular to that of the incident beam is 
detected; this is of particular importance in backscatter- 
ing, where we thereby eliminate the shortest scattered 
paths and favour the detection of long, deeply probing 
paths. The temporal evolutions of the speckle pattern are 
recorded at a frequency of 50 frames/s and 25 frames/s in 
transmission and backscattering, respectively, where we 
adjust the camera exposure time between 0.1-0.3 ms de- 
pending on the average intensity level. The speckle pat- 
terns are subsequently processed to obtain the temporal 
evolution of both the spatially averaged signal cj (t w ,r) 
from the transmission experiments and the spatially re- 
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(t w , t) that we present in dynamical 
While ci(t w ,r) is processed at different 
lag-times, we analyze cy (t w ,r) for the shortest accessi- 
ble experimental lag-time, r = 0.04 s, where we use the 
correction scheme outlined in section IV C of ref. [l8| to 
correct for the measurement noise. The size of the region 
of interest over which we take the pixel average is chosen 
to be 0.1 x 0.1mm 2 , which corresponds approximately to 
a bubble diameter at both ages. 
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FIG. 3: Panel displaying dynamical characteristics obtained 
from the far field transmission experiments at age 1 (left) and 
age 2 (right), (a) and (b): time resolved degree of correlation 
as a function of t w calculated for different lag times (a) from 
top to bottom: r = 0.04,0.2,0.68,1.6,4,12 s and (b) from 
top to bottom: r = 0.04,0.12,0.32,0.8,1.6,4 s. The solid 
and dotted vertical lines denote the moments at which the 
dynamical activity maps shown in Fig. [2] were acquired; the 
solid lines correspond to Fig. [2ja) and (c), respectively and 
the dotted line to Fig. [2jb) and (d). (c) and (d): time aver- 
aged correlation function as a function of lag time; to account 
for changes in the experimental geometry r is normalised by 
(L/l*) 2 . The continuous lines present the results of exponen- 
tial fits to the data, (e) and (f): probability density function 
of the 'instantaneous' decay time to (t w ) obtained by fitting 
the age dependent two-time correlation function to eqn (5). 
The continuous red lines present the sum of two Gaussian dis- 
tributions (green dashed lines) that are used to fit the data. 



III. RESULTS AND DISCUSSION 

Typical examples of the temporal evolution of the in- 
stantaneous degree of correlation obtained in our trans- 
mission far field experiments are shown for age 1 in 
Fig. [3Ja) and for age 2 in Fig. ^[b). The instantaneous 
degree of correlation systematically decreases with in- 
creasing lag-time, indicating that at any moment in time 
the system exhibits dynamical activity. Indeed, the in- 
stantaneous two-time correlation function can be read as 
a vertical cut through the data shown in Fig. [HJa) and 
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(b); the sporadic decreases in cj (t w ,r) observed at both 
ages thus effectively indicate that the dynamical activ- 
ity of our foams intermittently increases. Despite this 
variation, the overall time averaged correlation function 
cj (t w , r) exhibits a simple exponential decay, as shown 
in Fig. |31Jc) and l^d) for age 1 and age 2, respectively 
Such exponential decay is consistent with previous ex- 
periments 0, HI, [l9| and in principle would be consistent 
with the model described above, where the random inter- 
mittent restructuring of volume elements causes a decay 
of the correlation function that is indistinguishable to 
the one of continuously diffusing particles. However, as 
we will show, the decay of the mean correlation func- 
tion is determined by at least two dynamical processes, 
such that this model does not fully apply for the descrip- 
tion of the time averaged correlation function. Indeed, 
the inspection of the temporal evolution in the dynami- 
cal activity maps, shown as movie 1 (age 1) and movie 
2 (age 2) in the ESI [29[ reveals that the events mostly 
occur one at a time in the observational volume of our 
experiment. Moreover, the movies clearly reveal that the 
local restructuring of a volume element is a rather slow 
process. Indeed, a characterization of the event duration 
by calculating the temporal correlation of the dynami- 
cal activity at a given location reveals that the events 
typically lasts ~ 1 s at age 1 and ~ 1.8 s at age 2. 

Between events, the dynamical activity maps reveal 
that cf (t w ,r) is relatively homogeneous in space, as 
shown for both ages in Fig.^b) and Fig.^d). The direct 
comparison of the results obtained in our photon correla- 
tion imaging experiment to those obtained in the simulta- 
neously performed TRC-experiments shows that the fluc- 
tuations in the TRC-traces are effectively caused by sin- 
gle events of finite duration. However, during the event- 
free periods the spatially averaged dynamics probed in 
our transmission experiments is not suppressed. This 
can be seen for instance at the moments that are marked 
by a dotted line in Fig. [3ja) and Fig.[3][b), which are the 
moments at which the dynamical activity maps shown in 
Fig. [2ja) and (c) were acquired. Clearly, a (t w ,r) still 
fully decays during these event-free periods, which indi- 
cates that another dynamical process contributes to the 
dephasing of the scattered light. As this process does 
not lead to significant fluctuations in the degree of cor- 
relation, neither in space nor in time, we assume that 
this dynamics is quasi-continuous and occurs everywhere 
within the sample. Our findings are consistent with the 
findings obtained in recent speckle visibility experiments, 
where both the finite duration of the event and the exis- 
tence of an addition dynamical process were detected as 
well [l5|. However, in contrast to the speckle visibility 
technique, the TRC-technique allows us to assess the full 
lag-time dependence of the different dynamical processes 
contributing to the dephasing of the scattered light. As 
both the duration of the events and the event-free pe- 
riods between events are rather long compared to the 
experimental decay time of the time-averaged intensity 
correlation function, we can attempt to isolate the con- 



tribution of the quasi-continuous process by fitting the 
two-time correlation function with 

ci{tw > T) = ih cxp {-(^k)Y ] (5) 

where we allow the intercept l//3(t w ), the decay time 
To {t w ) and the stretching exponent p (t w ) to vary with t w . 
At both ages we find that the probability density func- 
tions (PDF) of the instantaneous decay times To (t w ) ex- 
hibit the same basic features: at shorter decay times, the 
PDF(t (t w )) exhibits a reasonably well-defined maxi- 
mum, at larger decay times the PDF(tq (t w )) is char- 
acterized by a shoulder, as shown for age 1 in Fig. [3je) 
and for age 2 in Fig. EJf), where we normalize the de- 
cay times with (L/l*) 2 to account for the varying ex- 
perimental conditions. Our results clearly show that the 
decay of the mean-correlation function is the result of 
two processes with different characteristic times. Indeed, 
we can attribute the time scale of the second maximum 
to the one characterizing the quasi-continuous process, 
while the time scale of the first maximum still consti- 
tutes a convolution of both the restructuring events and 
the quasi-continuous process. In order to determine the 
characteristic time describing the quasicontinuous pro- 
cess r c , we fit the PDFs to a pair of Gaussian distribu- 
tions that we find to describe the data reasonably well, 
as shown by the red continuous line in Fig.[3Je) and (f). 
From the resulting fit parameter To. maX 2(L/l*) 2 , we de- 
termine t c = 2 • TQ tmaX 2(L/l*) 2 and report the results 
in Table [TJ for better comparison to the original model 
outlined above (see eqn. and ©) we choose to con- 
sistently state decay times that apply to the decay of the 
field correlation function rather than the intensity cor- 
relation, which just entails to apply a factor of two to 
the decay time of the intensity correlation function. As 
compared to the mean decay-time r m = 2 ■ tq(L/1*) 2 de- 
termined from the exponential fit to the data shown in 
Fig. [3Jc) and (d), r c is larger by a factor of 1.2 and 1.5 
for age 1 and age 2 respectively. Interestingly, such mod- 
erate change in the decay time of the correlation function 
was also reported for experiments probing the foam dy- 
namics before and after a large oscillatory shear strain 
was applied [1] . Presuming that a large oscillatory shear 
strain erases internal stress imbalances, such that the re- 
arrangement events are temporally suppressed, it is likely 
that the residual dynamics observed after the application 
of shear corresponds to the quasi-continuous dynamics 
observed here between events. This would suggest that 
the quasi-continuous process is not directly linked to the 
release of imbalanced stresses. 

Finally, we note that the age dependence of the quasi- 
continuous process does not correspond to the one of 
the mean dynamics, as denoted by comparing the ratios 
T m (age 2)/r m (age 1) and r c (age 2)/r c (age 1) listed in Ta- 
ble [J This indicates that the event dynamics contributes 
to the mean dynamics so that the age-dependence is al- 
tered. Moreover, an inspection of the correlation be- 
tween tq (t w ) and p(t w ) reveals that the quasicontinu- 
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TABLE I: Comparison of coarsening and dynamical charac- 
teristics of Gillette foam at two different ages. t w age of the 
sample at the middle of the experiment; (at w ) mean bubble 
radius; d (at w ) /dt rate of bubble growth; dj/dt strain rate set 
by the coarsening process; r m characteristic decay time of the 
time and spatially averaged field auto-correlation function; 
r c decay time characterizing the quasi-continuous process; R 
rate of events per unit volume obtained from counting proce- 
dure; IV restructuring rate obtained from counting procedure; 
Td duration of intermittent events. 
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ous process is best described by a slightly compressed 
exponential, while during the occurrence of events the 
decay of the two time correlation function is described 
by a slightly stretched exponential. It is indeed remark- 
able that this combination between fast and stretched 
decays and slow and compressed decays leads to a sim- 
ple exponential decay in the mean correlation function. 
However, this high degree of convolution of event dynam- 
ics and quasi-continuous dynamics determining the mean 
dynamics significantly hinders the characterization of the 
rate of topological rearrangements and consequently their 
dependence on the coarsening process. 

By contrast, the characterization of the event dynamics 
is rather straightforward using the results of our photon 
correlation imaging experiments. Indeed, the temporal 
evolution of the events can simply be followed in movies 
like the one shown in the ESI [231 . J us * counting the 
events allows us to determine the event rate per unit vol- 
ume R reported in Table HI where we estimate the size 
of our observational volume as the product of the size 
of the field of view of our camera and the thickness of 
the cell. Clearly, the event rate within a given volume is 
significantly larger at age 1 than at age 2. However, to 
evaluate r r , the rate at which a volume element is effec- 
tively restructured, we need to account for the number 
of possible rearrangement zones within our observation 
volume, which scales with the age-dependent size of the 
rearranged zone £. An analysis of the event size, by calcu- 
lating the spatial correlation of cf^ (t w , r) introduced as 
d?4 (As, t) in ref. [lU, reveals that £ scales with the bub- 
ble size and extends over 10 bubble diameters, in agree- 
ment with previous assessments 

n m m. Th ° iar g° 

spatial extent and long durations of the events observed 



in our experiments are consistent with the events probed 
at the foam surface, where we observe that a series of 
neighbour-switching events, Tl-events [l[, rearranges a 
zone extending over ~ 4 bubble diameters. The even 
larger size of the event observed in diffusing wave spec- 
troscopy can be understood as that the event typically 
contains both a central volume element, in which the 
bubbles topologically reconfigure, and a halo of bubbles 
that respond elastically without rearranging, such that 
the zone that effectively contributes to the dephasing of 
light extends over ~ 10 bubble diameters @, laEl- 

To evaluate the restructuring rate, we thus assume that 
the size of the effectively topological rearranged zone is 
£ = 4 (a* J- As shown in Table [fl r r = R • |tt (4 (a t J) 3 
decreases by approximately a factor of 2.1 from age 1 
to age 2. To assess the relation between the topologi- 
cal rearrangements and the coarsening process, we eval- 
uate the strain rate set by the coarsening process as 
d-y/dt (AV/V) /At, with V and AV/At respectively the 
bubble volume and the change of the bubble volume per 
unit time at a given foam age. As shown in Table HI 
we find that the ratio of the age dependent strain rates 
[(d-y / dt(&ge I)) /(dj/dt (age 2))= 2.02] correspond, within 
the experimental error, to the ratio of the restructuring 
rates [T r (age l)/r r (age 2) = 2.1]. This strongly sug- 
gests that the events are a direct result of the coarsening 
induced strain between events. Estimating this critical 
strain as jy — (dj /dt)/T r , we find that j y sa 0.16 ± 0.01 
in reasonable agreement with the critical strain mea- 
sured for our foams in experiments probing the macro- 
scopic mechanical properties. This correlation between 
strain rate and event rate is not evident when consid- 
ering the decay rates of the mean correlation function 
t" 1 , which, following eqn. ([1]) [t^ 1 = R ■ |7r£ 3 ] , should 
effectively correspond to the restructuring rates if the 
dynamic light scattering signal would be determined by 
the random restructuring of volume elements only. Here 
the ratio r~ 1 (age l)/r~ 1 (age 2) would appear to cor- 
relate to the ratio of the bubble growth rather than to 
that of the strain rates, as shown in Table HI Indeed, 
the convolution of event-frequency, event duration and 
quasi-continuous dynamics that all exhibit a different 
age-dependence leads to an agc-dcpcndcnce in the mean- 
dynamics that is unsuitable for establishing relations be- 
tween the coarsening process and topological rearrange- 
ments. 

To address the question whether not the disappearance 
of bubbles, T2-events [l|, rather than the mean strain 
set by the coarsening process is at the origin of the me- 
chanical instabilities that lead to the intermittent rear- 
rangements observed, we estimate the number of bubbles 
disappearing within a time interval At as 

N b (t w )-N b (t w + A*) = -M - / o - — -3 ) 

(6) 

where ATj (t w ) and N b (t w + At) are respectively the num- 
ber of bubbles per unit volume at an age of t w and 
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FIG. 4: (a) Probability density function of time delay be- 
tween successive events. The solid line denotes the exponen- 
tial behaviour expected for the random occurrence of inde- 
pendent events, (b) Experimentally determined probability 
to find the (i + n) th event observed within our experimen- 
tal volume at a distance r from the i th event normalized 
by the probability expected for randomly distributed events, 
g(r) = PDF(r ) 

experiment 

/PDF(r) r 

andom ■ As g(r) varies 
smoothly with increasing n, g(r) obtained for n — 1-5 (filled 
squares), n = 6-10 (open squares), n = 11-15 (filled circles), 
n — 16-20 (open circles), and n = 36-40 (triangles) are aver- 
aged to reduce the statistical fluctuations. Inset: Probability 
that the (i + n) th event observed within our observational 
volume is located within the topologically rearranged zone of 
the i th event divided by the expectation value for randomly 
distributed events. 



t w + At. This estimate yields a rate for T2-events per 
unit volume that exceeds the rate of events per unit vol- 
ume R by a factor of = 10, such that it appears unlikely 
that T2-events are at the actual origin of the events ob- 
served. Finally, we note that single T2-events are small 
and involve the motion of only a few plateau-borders, 
which have been previously identified to be the princi- 
pal scattering sites in dry foams [3(| . T2-events are thus 
unlikely to be resolved in our experiments, as the resolu- 
tion is limited by the number of paths rearranged by an 
event divided by the total number of paths contributing 
to Cj (t w , t) of a given ROI. 

To further progress in our understanding of the topo- 
logical rearrangements, we characterize the temporal and 
spatial distribution of the events exploiting the tempo- 
rally resolved dynamical activity maps. We assess the 
time at which an event occurs as the moment where 
the dynamical activity is maximal during an event. At 
this moment, the location of the event is determined as 
the centre of the high dynamical activity zones, as those 
shown in Fig. 2(a) and (c). The results obtained from our 
analysis of the experiment at age 1 are shown in Fig. [4j 
qualitatively similar results are obtained for the experi- 
ment at age 2, but the quality of the data is somewhat 
lower as the number of events probed over the duration 
of our experiment is here lower than at age 1 . In agree- 
ment with previous results fl5| , we find that the time de- 
lay between successive events At e follows an exponential 
distribution, as shown in Fig. HJa). This is consistent 



with the Poissonian statistics expected for the random 
occurrence of independent events. To characterize the 
distribution of the events in space we determine the dis- 
tance, r, between the i th and (i + n) th event and calculate 
the probability density function of r, PDF(r) experiment, 
for a given n. To account for the limited geometry of our 
experiment we determine the PDF expected for a random 
distribution of events in space PDF(r) ranc i om using sim- 
ulations, where we implement the experimental bound- 
ary conditions. By normalising PDF(r) exper i ment with 
PDF(r) ran dom, we obtain the radial distribution func- 
tion g(r) = PDF(r) experiment / PDF(r) random , where we 
expect g(r) = 1 at all r for a random distribution of 
events in space. Instead, we find that g(r) > 1 for small 
r and n < 30, as shown in Fig.^Jb); as the decay of g(r) 
varies only slowly with increasing n we have here aver- 
aged each g(r) over a few values of n to reduce the sta- 
tistical fluctuations. The increased probability at small 
r indicates that the foam maintains a spatial "memory" 
of its dynamics; the occurrence of an event at a given 
location is thus correlated to previous rearrangements 
although no temporal correlation is observed. The inter- 
cept of g(r) decreases with increasing n, indicating that 
the structural memory effect weakens with time. More- 
over, we note that g(r) decreases to the background level 
at r corresponding to 4-5 bubble diameters; this indicates 
that the actual centre of an event is likely to be located 
within the topologically rearranged zone of a previous 
event. To quantify the time over which events are more 
likely to occur within such zone, we integrate g(r) for r up 
to 4 bubble diameters. The integrated value, p(0), com- 
pares the experimentally found event probability within 
an already rearranged zone of size r = 4 bubble diame- 
ters to that of randomly distributed events. It decreases 
slowly with n, reaching unity at n = 30, as shown in the 
inset of Fig.[4jb), where we have used a five-point moving 
average to smooth the data. Corresponding to n = 30, 
the spatial memory time of an event would thus extend 
up to w 30 s; using this time to evaluate the mean strain 
set by the coarsening process we find a strain that is 2 
orders of magnitudes lower than j y . 

These findings suggest that a fragile bubble configu- 
ration is likely to still be fragile after reconfiguration. 
Indeed, if we consider that the coarsening process leads 
to imbalances in the stress-configuration, an event will 
be set off where the excess stress exceeds the local yield 
stress. The reconfiguration process will proceed until the 
excess stress is again below the yield stress; this, how- 
ever, implies that the excess stress of the recently rear- 
ranged zone is just below the yield stress. As the coars- 
ening process continuously perturbs the stress configura- 
tions throughout the system, the likelihood that the yield 
stress within this zone is soon going to be exceeded again 
is then high, consistent with the observed behaviour. An 
increased probability to find an event within the already 
rearranged zone is also enclosed in the fact that the topo- 
logical rearrangement leads to an elastic deformation of 
the surroundings. Indeed, this elastic response is likely 
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to further contribute to the unfavourable stress balance 
within the just rearranged zone. These conditions de- 
note that coarsening induced rearrangements are fairly 
inefficient in relaxing unbalanced stress configurations 
and thus suppressing rearrangements. This contrasts 
with the effect of imposing large macroscopic oscillatory 
strains, which has been shown to efficiently suppress the 
occurrence of rearrangement events @, Indeed, im- 
posing a strain that significantly exceeds the yield strain 
leads to a reconfiguration of the ensemble of the bubbles; 
both the initiation and cessation of bubble rearrange- 
ments are here not governed by threshold configurations; 
moreover, the topological reconfiguration of the ensemble 
eliminates the effect of the elastic deformation, such that 
fragility is not preserved. 

Finally, we note that a persistence of stress-imbalances 
within a given region has been recently inferred from ex- 
periments probing the spatial heterogeneities in the dy- 
namics of compressed microgel systems [1[ . There it was 
observed that well-localized zones of high dynamical ac- 
tivity coexisted with zones of low dynamical activity and 
that the high dynamical activity within a zone persisted 
in time, which indicates that rearrangements not neces- 
sarily lead to a balanced stress configuration. 

To some extent this behaviour is reminiscent of dy- 
namical facilitation in glasses, where an immobile zone 
can only become mobile when the neighbouring zone is 
already mobile [3l| • As a direct consequence of the highly 
cooperative nature in the dynamics of glassy systems the 
onset of mobility in a given zone correlates here to the 
mobility of the neighbouring zone. In foams, the onset 
of mobility in a given zone correlates to previous mobil- 
ity within the same zone. One may speculate that both 
correlations are related to the actual origin of the dy- 
namics in glasses and deeply jammed systems. While the 
thermally-driven dynamics in glasses imposes cooperativ- 
ity, i.e. fluidization of a neighbouring region, rearrange- 
ments driven by unbalanced stresses in deeply jammed 
systems increase the 'degree of jamming' in the surround- 
ings, thereby inducing an imbalance of internal stresses 
at the location that has just been rearranged. 

IV. CONCLUSIONS 

In conclusion, we have investigated spatial and tempo- 
ral heterogeneities in the dynamics of a coarsening foam 
at two different ages, conducting simultaneously experi- 
ments in the transmission and backscattering geometry 
of diffusing wave spectroscopy. In transmission, we per- 
formed the experiments in the traditional far field limit, 
while we used the photon correlation imaging configura- 
tion [l6j in the backscattering experiments. The com- 
bination of these experiments allows us to critically ad- 
dress different contributions to the dynamic light scatter- 
ing signal. In agreement with previous experiments [l5j . 



we find that the foam dynamics is determined by two 
processes: intermittent bubble rearrangements of finite 
duration and a quasi-continuous process that is spatially 
homogeneous. Both processes contribute to the dynamic 
light scattering signal, which significantly complicates 
the interpretation of the mean signal. For the charac- 
terization of the event dynamics, it becomes evident that 
the photon correlation imaging technique is a particu- 
larly useful tool, as it provides the means to assess both 
the frequency and the position of the events without the 
need of modelling the actual magnitude of the dynamic 
light scattering signal. These experiments reveal that 
the event frequency is directly correlated to the strain- 
rate imposed by the coarsening process, which indicates 
that the events are triggered by a well defined strain. 
Our analysis of the spatial correlation between succes- 
sive events uncovers that the events are not randomly 
distributed in space. Indeed, we find that the probabil- 
ity to find an event at a location that has been previously 
rearranged is significantly higher than expected for a ran- 
dom spatial distribution of events. This unexpected be- 
haviour indicates that a locally unstable zone rearranges 
to a state that is just below the yield stress, such that 
the likelihood that this zone is soon going to be again 
rearranged is high. Such persistence might be a general 
characteristic of stress-driven dynamics that could be ob- 
served in other jammed systems where the direct contact 
between the constituents may result in an unbalanced 
stress-configuration. 

Our findings evoke several interesting questions that 
could be addressed in future experiments. In particular, 
the experiments probing the effect of a macroscopic ap- 
plied strain could be used to explore the fragility of a 
rearranged zone. One could envision that a small oscil- 
latory strain would facilitate the secondary events, such 
that the events would merge within the typical duration 
of a single event; in such a case it is likely that the spatial 
distribution of events would appear to be random. More- 
over, it would be interesting to investigate the evolution 
of the spatial heterogeneities in the dynamics of a foam 
that has been subjected to a large external shear, as this 
has been done for the mean dynamics in ref . @, • Such 
investigations should help to assess the efficiency of re- 
structuring events in equilibrating stress imbalances and 
could lead to a better understanding of the quasicontin- 
uous process whose origin remains unclear. 
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